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ABSTRACT 

We study the UV spectra of narrow-line Seyfert 1 (NLSyl) galaxies and 
compare them with "normal" AGN. Similar to their optical lines, the NLSyls 
show narrower UV lines. They are also characterized by weaker CIVA1549, 
CIII] A1909, and stronger A1IIIA1857 emission. These UV line properties add 
to the optical and X-ray properties known to be part of the Boroson & Green 
eigenvector 1. We show that the steep soft-X-rays, which characterize the 
NLSyls SEDs, change the equilibrium of the two phase cloud-intercloud medium 
resulting in somewhat higher BLR cloud densities, lower ionization parameter 
and larger BLR radii. These modified conditions can explain the unusual 
emission line properties we find in NLSyl. 

Using a specific model of an accretion disk with corona presented by Witt, 
Czerny & Zycki, we also show that the steep soft and hard- X-ray continua can 
be explained if the L/LEdd ratios are larger than in "normal" Seyfert 1/Q SO 
strengthening earlier suggestions that the L/L Edd is the physical parameter 
driving this eigenvector. 



Subject headings: galaxies: active — galaxies: Seyfert — accretion, accretion 
disks 
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1. Introduction 

Narrow-line Seyfert 1 galaxies (NLSyl), first suggested as a distinct class of AGN by 
Osterbrock and Pogge (1985), are characterized by Balmer lines whose FWHM is smaller 
than typical Seyfert 1 galaxies i.e. 500 < FWHM < 2000 km s -1 , slightly broader than 
the forbidden lines. On the other hand they are clearly different from Seyfert 2 galaxies 
since the ratio of [OIII] A5007 to H/3 is <3, i.e. below the limiting value found by Shuder 
and Osterbrock (1981) to discriminate between Seyfert 1 and Seyfert 2 galaxies. In NLSyl 
strong Fe II optical multiplets and higher ionization iron lines (e.g.[FeVII]6087A and [FeX] 
6375A) are often present. These are usually seen in Seyfert 1 and not in Seyfert 2 galaxies. 

Many NLSyls have an unusually strong big blue bump (BBB) which, when compared 
to typical Seyfert 1 and QSO BBBs, is shifted towards higher energies, sometimes even 
out of the optical/UV range (at least one object actually peaks in the soft X-ray band: 
RE J1034+396, Puchnarewicz et al. 1995). Its high frequency tail is clearly seen in soft 
X-rays and these objects have generally steeper soft X-ray continua than is "typical" 
for Seyfert 1 galaxies (Boiler, Brandt & Fink 1996), meaning that they have a stronger 
soft-X-ray excess over the hard X-ray power law. The intrinsic hard-X-ray continua of 
NLSyls are also generally steeper (Brandt, Mathur, Elvis 1997) than in typical Seyfert Is. 
The NLSyls are usually only weakly absorbed in the soft X-rays (Boiler, Brandt & Fink 
1996) and in many cases both the UV flux and soft X-ray flux are strongly variable. 

NLSyl objects are generally radio-quiet and their radio powers are typical of those 
found in other Seyfert galaxies (Ulvestad, Antonucci & Goodrich 1995). 

There is no widely adopted view on the basic reason why the continua of NLSyl 
galaxies are different from classical Syl. The two most probable explanations of the stronger 
big blue bumps in these objects are pole-on orientation (Puchnarewicz, Mason, Cordova 
1994, Wilkes 1998), and higher accretion rate relative to the mass of the central object 
(e.g. Boiler, Brandt, Fink 1996; Wandel 1997, Czerny, Witt & Zycki 1997, Pounds, Done 
& Osborne 1995). Steeper hard X-ray spectra and strong permitted Fell lines are possibly 
a secondary effect of the atypical shape of the soft-X-ray continuum (Pounds, Done & 
Osborne 1995, Brandt, Mathur & Elvis 1997; Wilkes, Elvis and McHardy 1987, Shastri et 
al 1993). 

Wilkes et al. (1999), studying a sample of low redshift quasars and Syls, and their 
relations between optical/UV emission lines and the continuum, found that the four NLSyl 
in their sample show smaller equivalent widths of CIII] and CIV lines than typical AGN 
(EW(CIV) < 40A for NLSyl, while 30A < EW(CIV) < 200A for other AGN). In this paper 
we investigate in detail the UV line properties of a sample of NLSyl objects to determine 
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whether the weakness of the carbon lines is typical of these objects, constituting an 
additional property which distinguishes them from "normal" Syl galaxies. We investigate 
the physical conditions of the BLR which may explain these systematic differences. We 
also discuss the possibility that these objects have luminosities close to their Eddington 
luminosity. 

2. UV line measurements 
2.1. The sample 

From the currently known set of NLSyl (Boiler, Brandt & Fink 1996, Greiner et al. 
1996, Puchnarewicz et al. 1992, Puchnarewicz et al. 1994, Brandt, Fabian & Pounds 1996, 
Grupe et al. 1996 , Moran et al. 1996, Brandt - private communication, Wilkes et al. 1999) 
we have defined a subset of 11 objects (Table 1) for which UV spectra are available either 
from the HST (5 objects) or IUE archives. The IUE spectra were taken from Lanzetta et 
al. (1993), and the reduced HST spectra from Dobrzycki (private communication, see also 
Bechtold et al. 2000). 

2.2. Line parameters 

We have measured the EW (Table 2), line ratios (Table 3) and line widths (FWHM; 
Table 4) of all prominent UV lines: LyaA1216, CIVA1549, CIII] A1909, Silll] A1892, 
A1IIIA1857, SilV+OIV] A1400 blend, and MgIIA2798. 

The EW and FWHM of IUE spectra were measured using the splot task in IRAF: the 
EW by fitting a linear continuum to the data and integrating across the observed emission 
line (keystroke 'e'), the FWHM by measuring the width at half the flux in the line peak 
above the continuum. The same procedure was applied when the line parameters in the HST 
spectra were measured, although a different program findsl (provided by Aldcroft, Bechtold 
& Elvis 1994), specially written to handle the HST data, was used. The line parameters 
presented in Tables 2,3,4 have been corrected for absorption: for weak absorption by using 
a linear fit across the absorption line, for strong absorption (as in PG 1351+560 and 
PG 1411+442) by assuming a symmetric emission line profile and reflecting the unabsorbed 
wing about the peak. 

As has been noted by Vestergaard & Wilkes (2000), CIII] is blended with Felll UV34 
A1914 line, which should be taken into account especially when the CIII] line is weak. In 
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three spectra (IZwl, PG 1211+143, Mrk 478), where the Felll UV34 was clearly visible, 
we subtracted this line (modeled as a Gaussian centered at A=1914A rest frame) from the 
CIII] blend. 

2.3. Comparison of NLSyl with "normal" AGN 

In this section we compare the UV line properties of NLSyl with Seyfert 1 galaxies 
and quasars. Fig. la shows the EW of Lya, CIV and Mgll of our NLSyl sample (shaded 
areas) compared to the sample of Seyfert 1 galaxies (dotted line) from Wu at al. (1983) 
(their sample includes three NLSyl: IZwl, Mrk 478, IIZwl36 which we excluded here) and 
low redshift quasars from Wilkes et al. (1999), Corbin & Boroson (1996), and radio-loud 
quasars from Baldwin, Wampler & Gaskell (1989), combined and denoted by a dashed line. 
The EW of CIV and Mgll lines are significantly smaller in NLSyl than in the broad line 
Seyfert 1 galaxies and quasars. The K-S test yielded a 0.001 chance that the EW of CIV 
and Mgll in NLSyls and Syls are drawn from the same population. For Lya the chance 
was < 0.02. When compared to QSO the significance remained strong for CIV (p < 0.01) 
and Mgll (p < 0.025), while for Lya the distributions are similar (p > 0.5). The smaller 
EW of the carbon and Mgll lines cannot be due to a simple continuum increase, as this 
would effect the EW of all lines equally, while EW(Lya) is not significantly smaller. 

As can be seen from Table 3 the CIV/Lya (mean 0.25±0.09), CIII]/Lyct (mean 
0.05±0.05) and the Mgll/Lya (mean 0.05±0.03) ratios are smaller compared to those 
typically observed in Seyfert 1 galaxies (Wu et al. 1983 give observed ranges: 0.35-2.01, 
0.03-0.39, 0.07-0.63 respectively) and quasars (observed range: 0.3-1.04, 0.15-0.3, 0.15-0.35). 
To show this more clearly Figure 2 shows the CIII]/Ly« vs CIV/Lya line ratios for NLSyl 
in our sample (denoted as filled squares) with the Seyfert 1 sample (denoted as circles) and 
quasars from Laor et al. (1995), Christiani & Vio (1990), Wilkes et al. (1999) and narrow 
line quasars from Baldwin et al. (1988). 

As the lines in NLSyl are narrow, we can clearly resolve the components of the 
CIII] +SiIII]+ A1III blend (especially in the HST data, in the IUE data the S/N is often too 
low). From Table 3 it is also clear that the Silll] line in most of the NLSyl is very strong 
compared to the CIII] line. Also the SilV+OIV] blend is strong compared to CIV (mean 
SiIV+OIV]/CIV ratio in NLSyl is 0.49±0.26, larger than the mean ratio of 0.3 in quasars 
from Francis et al. 1991). However, the SiIV+OIV]/Lya; ratio is in the range of normal 
AGN, indicating that the large SiIV+OIV]/CIV ratio is due to weaker CIV emission. The 
A1III doublet in NLSyl is rather strong (equivalent width ~ few A - see Table 2). 
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Fig. 1. — a) Comparison of Lya, CIV, and Mgll equivalent widths in our NLSyl sample 
(shaded areas) with the Seyfert 1 sample from Wu et al. (1993) (dotted line) and QSOs 
(dashed line) from Wilkes et al. (1998), Corbin & Boroson (1996), and Baldwin, Wampler 
& Gaskell (1989) 

Although broader than H/5 the UV lines in NLSyl are narrow compared to other AGN. 
In Table 4 and Fig. lb we compare our objects with samples of low redshift quasars from 
Corbin & Boroson (1996) and Wilkes et al. (1999), and with a radio-loud sample from 
Baldwin, Wampler & Gaskell (1989). 

3. Discussion 

We will now investigate what the line strengths and line ratios tell us about the physical 
properties in the BLR clouds of NLSyl (Section 3.1). Then we will discuss the continuum 
properties of NLSyl (Section 3.2) and investigate what they indicate about their central 
engine. In conclusion we show how the deduced differences between the central engines of 
NLSyl and "normal" AGN can explain their different, observed emission line spectra. 
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Fig. 1. — b). Comparison of FWHM of Lya, CIV, and Mgll full widths at half maximum, 
with the same coding as in Fig la. 



3.1. Physical properties of the BLR clouds in NLSyl 

More than ten years ago Gaskell (1985) noticed that Seyfert 1 galaxies with narrow 
H/3 lines of FWHM < 1600 km s -1 show lower H/3 equivalent widths than typical Seyfert 1. 
He interpreted this finding as a result of collisional destruction of H/3 in the higher density 
BLR clouds in these objects. Although we do not study the optical spectra of NLSyl in 
this paper, we will now investigate whether the UV spectra lead to a similar conclusion. 

Rees et al. (1989) calculated the line intensities for different BLR cloud densities at 
constant column density (10 23 cm -2 ) and ionization parameter (U = 10~ 2 ). They found 
that optically thick lines such as hydrogen and carbon lines have a fairly constant intensity 
up to a certain density, above which these lines become thermalized and their intensity 
drops considerably. For hydrogen lines, CIV and Mgll this critical density is ~ 10 10 cm~ 3 . 
For the semi-forbidden lines CIII] and Silll] it is around 5 x 10 9 cm~ 3 and 10 11 cm -3 
respectively. The Lya and CIV lines are usually strong coolants at densities smaller than 
these critical values, but as the density increases and these lines become thermalized, other 
high-excitation lines such as CIIIA977 and A1IIIA1857 take over the cooling. 



In the previous section we showed that the UV spectra of NLSyls, when compared to 
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Fig. 2. — Comparison of CIIIj/Lyo; and CIV/Lya ratios for NLSyl objects (filled squares) 
and Seyfert 1 galaxies from Wu et al. (1983, circles) and QSOs from Laor et al. (1995, stars), 
Christiani & Vio (1990, crosses), Wilkes et al. (1998, triangles) and narrow line quasars from 
Baldwin et al. (1988, diamonds). (For CIII] we used the sum of CIII]+SiIII]+AlIII to allow 
comparison with other samples where the broader lines prevented the authors from separating 
these components.) 



"normal" Seyfert 1 and QSO galaxies, show weaker carbon and Mgll lines. Although the 
wavelength of our spectra does not cover the range of CIII A977, the A1III A1857 doublet is 
clearly seen (where the S/N is high enough) and is especially strong in IZwl (see Table 2). 
All these line properties suggest that in NLSyl objects the BLR clouds have higher densities 
than the BLR clouds in "normal" AGN. We will estimate how much higher by studying the 
line ratios in the following section. 



a) The line ratios 



The CIII] and CIV to Lya ratios are often used as a density indicator. This is because 
the carbon lines are collisionally excited (hence sensitive to density), while Lya is not (note 
however that Mathur et al. 1994 showed that for large ionizing parameters, U > 0.1, where 

U = , 1Ry 2 d hu , CIIIl ceases to be a density indicator). These line ratios are also a sensitive 



- 8- 



function of the ionization parameter U, therefore we investigate the relation of these line 
ratios to density and U. 

We have calculated line ratios using the photoionization code CLOUDY (version 80.07, 
for reference see Ferland 1991). First as an input ionizing continuum we took a standard 
AGN continuum (table agn; Mathews & Ferland 1987). With this continuum, the observed 
lines were only reproduced with higher densities in the CIII] than the CIV emitting clouds. 
This requires a steep increase in cloud density with radius, which is contrary to expectations 
and seems unrealistic. Then we used the spectral energy distribution (SED) of the NLSyl 
PG 1211+143 as the ionizing continuum. Although no detailed study of the SEDs of 
NLSyls has been made and is beyond the scope of this paper, PG 1211+143 is typical 
of those studied to date with a x = 2.13 ± 0.22, where F v = v~ a (Wang, Brinkmann & 
Bergeron 1996; where typical NLSyl slopes are in the range 1.5 to 3.5 - see Boiler, Brandt 
& Fink 1996) and ai = 0.90 (i.e. slope measured between 1 /im and 2500A, where typical 
NLSyl values are 0.4 to 2.8 - see Lawrence et al. 1997). The IR to hard X-ray SED of this 
object was taken from Elvis et al. (1994) and is reproduced here in Fig. 3. The SED was 
linearly interpolated between the observational points in the optical/UV region. The EUV 
continuum was determined by a linear interpolation between the lowest energy point in the 
X-ray range and the highest in the UV, providing a conservative (i.e. low) estimate of the 
number of EUV photons. We investigated a range of cloud densities (n(H) = 10 8 — 10 13 
cm -3 , well within the range of the applicability of the photoionization code CLOUDYQ) and 
ionization parameters (U = 10~ 3 — 10 , where U = 10~ 2 is the value for the "standard" 
BLR - Davidson & Netzer 1979). The metal abundances were assumed to be solar and the 
cloud column density 10 23 cm -2 . The calculated line ratios are plotted in Fig. 4, where 
the observed line ratios for our NLSyl are denoted by horizontal lines. The CIV/Lya 
ratio (Fig. 4a) depends very strongly (more than any other line ratio) on the value of the 
ionization parameter as well as on the density of the BLR clouds. Small values of the 
ionization parameter (U = 10~ 3 ) are clearly favored by our data, for which densities of the 
order of 10 11 cm -3 to 10 12 cm~ 3 are needed to produce the low observed CIV/Lya ratios. 
For the same, small ionization parameter the observed CIII]/Lya ratios indicate densities 
between > 10 9 cm -3 and 10 11 cm -3 (see Fig. 4b), which are smaller than the cloud densities 
inferred from the CIV/Lya ratios. This strongly suggests that the CIII] and CIV lines are 
formed in different clouds, which implies a stratified BLR. This is as expected from the 



3 "The hydrogen and atoms and ions of helium are treated in the code as 10-level atoms. The treatment of 
the heavy elements is not as complete as hydrogen and helium, but a 3-body recombination is included as a 
general recombination process. [...] The physical high-density limit is set by the approximate treatment of the 
three-body recombination-collisional ionization (< 10 13 cm~ 3 ) for the heavy elements and the approximate 
treatment of line transfer". - see Ferland (1991). 
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Fig. 3. — The infrared to hard X-ray spectral energy distribution of PG 1211+143 (from 
Elvis et al. 1994), which was used as an input ionizing continuum in the photoionization 
calculations. The dotted line indicates the standard AGN continuum of Mathews & Ferland 
(1987) used in CLOUDY normalized at 1/im. The NLSyls have a more pronounced BBB 
and a strong soft-X-ray excess. 

results of reverberation mapping (e.g. NGC 5548 Korista et al. 1995, NGC 7469 Wanders 
et al. 1997, Fairall 9 Rodrigues-Pascual et al. 1997) which show that the CIV and Lyct line 
fluxes vary with a smaller time delay, relative to the UV continuum, than the CIII] lines, 
indicating that CIV, Lya emitting clouds lie nearer to the central engine than the CIII] 
emitting clouds. The density of the CIV and Lya emitting line region is typically estimated 
(e.g. Peterson et al. 1985) to be ~ 10 11 cm~ 3 , while the CIII] region is ~ 10 9 ' 5 cm" 3 . Thus 
the density of the CIV, Lya emitting clouds in our NLSyls is comparable or somewhat (< 
10 times) larger than in normal AGN. 

To further constrain the density of the CIII] emitting clouds, we need to investigate 
its ratio to a line which is formed in the same clouds for example Silll]. For U = 10~ 3 the 
density of the emitting gas inferred from the Silll] / CIII] line ratio is between 10 9,5 cm -3 
and 10 10 - 5 cm~ 3 (see Fig. 4c, we omit here PG 1411+442, which is a BAL QSO and has 
n(H) ~ 10 9 cm" 3 ). The ratio of SiIII]/CIII] is larger than the typical value of ~ 0.3 ± 0.1 
seen in quasars (Laor et al. 1995) for all objects (except for BAL QSO PG 1411+442; see 
Table 3 and Fig. 4c) This high ratio is probably the result of the suppression of CIII] while 




9 10 11 12 

log n(H) 



Fig. 4.— Calculated line ratios of a) CIV/Lya, b) CIII]/Lya, c) Silll] /CIII] , d) 
SiIV+OIV]/Lyo;, e) Mgll/Lya, for several different ionization parameters. Triangles denote 
\ogU = —1, squares log £7 = —2, filled circles logU = —3 and stars logU = —2.5. Horizontal 
lines show the observed line ratios. 

Silll] remains strong, due to the smaller critical density for CIII] (> 5 x 10 9 cf. 10 11 cm" 3 
for Silll], see Section 3.1a). 

In Fig. 4d we present the observed and calculated SiIV+OIV]/Lyo; ratios. For U = 10~ 3 
the clouds have density of the order of 10 11 cm -3 to 10 12 cm" 3 , similar to the range for CIV 
and Lja emitting clouds. 

The Mgll/Lya ratio, on the other hand, is very small, and cannot be reproduced by 
clouds with small ionization parameter even when the column density, Nh, is varied over 
the range 10 22 - 10 24 cm" 2 . Only larger than standard ionization parameters, in the range 
1CT 1 to 1CT 2 , can produce the observed ratios (Fig. 4e). This may indicate that the Mgll 
lines do not form in the same region as the other BLR lines. We will return to this problem 
in the next section. 



To summarize the results of this subsection we conclude that the unusual UV line 
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ratios in the NLSyl objects can be explained if the BLR clouds have 10 times lower ionizing 
parameters (logC/ ~ —3) and a few times (<10) higher densities (n(H)~ 10 n ~ 12 cm -3 for 
Lyct, CIV, SilV emitting clouds and io 9 - 5-10 - 5 cm -3 for CIII], Silll] emitting clouds) than 
normal AGN. The BLR is clearly stratified with CIV, Lya, SilV producing clouds lying 
closer to the central engine, while CIII] and Silll] emitting clouds lying predominantly 
further out. The Mgll emission cannot be produced by the same cloud population, 
suggesting that these lines form in a different region. 

Although clouds with a wide range of properties are likely to exist in the broad-line 
region, it was shown by Baldwin et al. (1995) and Korista et al. (1996), that each emission 
line is most efficiently produced in gas with the optimum parameters for that line. These 
are the so called locally optimally emitting clouds or LOCs. Thus our modeling derives 
the parameters of the LOCs for each line so that the line fluxes and ratios provide a good 
approximation to a detailed multi-zone model of the BLR (see Baldwin et al 1995, Korista 
et al. 1996 and our Table 3), which is beyond the scope of this paper. 

b) The weak Mgll problem 

The Mgll line is surrounded on both sides by Fell emission. As the Fell emission in 
NLSyls is usually very strong (Boiler, Brandt & Fink 1996), it is possible that the wings of 
Mgll disappear in the stronger iron bumps. This effect could lead to an underestimation of 
the Mgll emission of up to a factor of two in the strongest Fell sources (such as IZwl - see 
Vestergaard and Wilkes 2000). However, even if our measurements were underestimating 
the Mgll emission by such a large factor, the real Mgll/Lyo; ratio would still be much 
smaller than that observed in quasars or in the lower end of the range for Seyfert Is (see 
Table 3). The ionization parameter U inferred from the observed line ratio would be 
~ 10 _1 — 10~ 2 (see Fig. 4e), still larger than that inferred from the other emission lines. 

Thus we conclude that the Mgll emitting clouds have a different value of the ionization 
parameter, and are formed in a physically different region of the BLR (consistent with the 
different time lags shown by Lja, CIV and Mgll lines in NGC 5548). Photoionization 
models predict that Mgll is either formed in a partially ionized zone (PIZ) of the BLR 
clouds or in a low ionization region (LIL) separate from the high ionization region (HIL) 
where the Lya, CIV, CIII] lines are formed ( Collin- Souffrin et al. 1988). If the Mgll line 
is formed in a PIZ it is possible that the stronger BBB in NLSyls will push the ionization 
front further back into the cloud, resulting in a smaller PIZ and weaker Mgll emission, than 
for objects with a "normal" BBB. 

However, the weaker Mgll emission is not consistent with the stronger Fell optical 
emission observed in NLSyl, which in photoionization models is predicted to be formed 
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in the same region (PIZ, Krolik & Kallman 1988, LIL, Collin-Souffrin et al. 1988). This 
inconsistency suggests that the Fell emission is instead generated in a different region from 
Mgll. The observations of line variability in NGC 5548 (Sergeev et al. 1997) showed that 
the Fell optical multipets have a very long time lag of several hundred days, while the Mgll 
has a 30-50 day time lag, also implying that these lines are formed in different regions. 
The Fell lines may be produced in the outer regions of the accretion disk as suggested 
by Dumont & Collin-Souffrin (1990) or in a separate, mechanically heated region closely 
related to the compact radio source as in Collin-Souffrin, Hameury & Jolly (1988) (hence 
the observed anti-correlation of Fell emission and radio flux). 

3.2. High luminosity to the Eddington luminosity ratio 

It has been suggested by a number of authors that NLSyl galaxies as a class have 
systematically higher ratios of their luminosity to the Eddington luminosity, i.e. they 
have systematically lower masses in a given luminosity range than Syl galaxies and QSOs 
(Pounds, Done & Osborne 1995, Wandel 1997). This suggestion was made based on the 
analogy with the Galactic black hole candidates. We will address this suggestion now. 

3.2.1. Continuum properties 

One of the current explanations of the soft-X-ray excess in AGN is reprocessing of the 
hard X-rays by partially ionized, optically thick matter, probably in the accretion disk. 
The model describes well the soft X-ray continuum of low-luminosity, flat a so f t x Seyfert 
galaxies, but has problems with fitting the steepest a so f t x spectra (see Fiore, Matt & 
Nicastro 1997), which characterize NLSyl. The steep a so f t x can instead be explained by 
emission from the innermost part of an accretion disk which is then Comptonized by an 
optically thin, hot corona surrounding the disk (Czerny & Elvis 1987; Laor et al. 1997). 

Theoretical models which can explain both the presence of the BBB and the hard X-ray 
emission are based either on radial or horizontal stratification between the hot optically thin 
and cold, optically thick accretion flow (Wandel & Urry 1991, Shapiro, Lightman & Eardley 
1976, for a review see Wandel & Liang 1991). In this paper we use the model of an accretion 
disk corona (ADC) by Witt, Czerny & Zycki (1997), where the corona itself accretes and 
generates energy through viscosity, and the division of the flow into optically thin and 
optically thick regions results from the cooling instability discussed by Krolik, McKee and 
Tarter (1981). Such a model is able to predict the fraction of the energy generated in the 
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corona instead of adopting this quantity as a free parameter. The model is fully defined 
by 3 parameters: the mass of the central black hole (M^), the accretion rate or the ratio 
of the luminosity to the Eddington luminosity {L/LEdd) and the viscosity parameter (a V i S , 
assumed to be the same in both the disk and the corona). The model predicts a systematic 
change in the opt/UV/X-ray spectral energy distribution due to a change in L/LEdd- A 
larger ratio results in a more pronounced BBB, which is shifted towards higher energies 
(resulting in stronger soft-X-ray emission and hence steeper soft X-ray slopes). 

We have determined continuum properties predicted by this model over a large range 
of L/LEdd (0.001 to 0.7), a vis (0.02 to 0.4) and black hole masses (10 6 to 1O 1O M ). We 
then compared the observed continua of our NLSyl with the UV luminosity at 2500A, and 
the soft and hard X-ray slopes (a so f t x and ahardx) predicted by the model. Table 5 shows 
the observed a so f t x (from ROSAT) and ahardx (from ASCA) slopes for each object, while 
Table 6 gives the best fitted model parameters for each object. Our model was able to 
reproduce the steep soft and hard X-ray slopes within the observed uncertainties for most 
of the NLSyl. However for two objects (IZwl, PKS 0558—504) we did not succeed in fitting 
both the soft and hard-X-ray slopes simultaneously. This may be due to the way we treat 
Comptonization in our model (see Janiuk & Czerny 1999 for further details). In Fig. 5a,b 
we show how the X-ray slopes change with the model parameters. Each curve represents 
one value of L/LEdd and a V i S and a full range of black hole masses, where smaller Mbh lie 
at smaller 2500A luminosities. We see clearly that only the large ratios of L/ L E dd can give 
the steep, observed soft X-ray slopes. 

As has been shown by Czerny, Witt & Zycki (1997), quasars radiate usually at ~ 
0.01-0.2 of their Eddington luminosity, while Seyfert galaxies radiate at ~ 0.001-0.3. Our 
NLSyl (where we use the same ADC model as Czerny, Witt & Zycki 1997 to fit the 
parameters of the central engine) radiate at L/L E dd ~ 0.27-0.58, much larger than the 
typical AGN. The masses of the central black hole calculated from the model (~ 1O 8 M to 
10 9 M Q ) for our objects are of the same order as masses found in typical Seyfert 1 galaxies, 
but the bolometric luminosities are larger, and comparable to those of QSOs (see Wilkes et 
al. 1999, Table 12 for comparison). This is deduced from the stronger, higher energy BBB 
and places the NLSyl in a transition zone between the Syls and QSOs i.e. among Syls with 
larger luminosities or QSOs with lower masses. We note that, while the absolute numbers 
we deduce depend upon the particular ADC model used, the general trends do not. 
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42 44 46 48 42 44 46 48 

log(i/L„[2500A]) log(i>L„[2500A]) 



Fig. 5. — Comparison of the a) soft and b) hard X-ray spectral slopes of NLSyl objects in our 
sample with those predicted by the ADC model. Each curve represents one value of Lj ' L E dd 
and a V i S with varying black hole masses as labeled (smaller lie at the smaller 2500A 
luminosities). Dashed lines denote small L/L Ed d=0-01, solid line L/L Ed d=0.3. The value 
of viscosity parameter are as labeled. Clearly the larger L/L Edd are needed to reproduce 
the steep soft X-ray slopes of NLSyls. The following symbols denote: open triangle - IZwl, 
filled circle - PKS 0558-504, filled triangle - PG 1211+143, filled square - IRAS 13349+248, 
filled diamond - PG 1351+640, open circle - Mrk 478, star - IIZwl36. 

3.2.2. Density and the radius of the BLR 

The structure and the dynamics of the BLR is complex, as suggested by variability 
studies in the case of Seyfert galaxies. However, we can analize the scaling properties of the 
whole BLR of an object with the properties of the central source, including the shape of the 
X-ray continuum. 

The BLR gas Compton heated by the ionizing continuum will form (in any geometry) 
two phases: a cool phase with T c ~ 10 4 (the BLR clouds) and a hot phase with T h ~ 10 8 
(the intercloud medium, see Krolik and Kallman 1988, Czerny & Dumont 1998, Wandel & 
Liang 1991), when in equilibrium. The precise values of these temperatures depend on the 
shape of the continuum. 
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In the context of the two-phase model, we will now investigate how the properties of 
the BLR change due to the steeper X-ray continuum of a NLSyl. We use the ionization 
parameter of Krolik, McKee & Tarter (1981): 

r-n 2.3i^j ori 2.3-F 1 j on . . 

" ~~ rr) ~ ckp c T c V / 

where p is the total pressure, p c and T c the density and temperature of the cold phase, and 
F ion is the flux above 1 Ryd determined by the ionizing luminosity of the central source 
L ion and the current radius r (where effects of geometry have been neglected): 

Fion = Uonl^r 2 . (2) 

The two phases coexist at a value of the ionization parameter, E^, which scales with 
the hot phase temperature, T h in the following way (Begelman et al. 1983): 

The BLR is most probably radially extended. For the purpose of exploring the various 
dependencies, we determine a representative radius for the BLR. Note that this is a scaling 
factor rather than the specific radius at which a particular emission line is generated. If the 
cloud number density profile is flatter than r~ 2 , then most of the emission would come from 
the outer radii of the BLR. As in the case of the Inverse Compton heated coronae discussed 
by Begelman et al. (1983), a nearly hydrostatic corona will exist up to a radius where 
the temperature of the hot medium is equal to the "escape" temperature (i.e. the virial 
temperature). At larger radii the corona is heated to temperatures exceeding the escape 
temperature, becomes unstable and forms an outflowing wind. We therefore identify the 
outer edge of the BLR, r BLR with the radius where the hot medium temperature is equal 
to the virial temperature 

GM bh m H 

kT h = (4) 

tblr 

The size of the BLR expressed in units of the Schwarzschild radius, Rschw is then given 

by: 

f'BLRl 'Rschw = n]rr > (5) 

2kT h 

so a lower value of the hot medium temperature in NLSyl galaxies is consistent with larger 
values of r BLR jR Schw and, consequently, lower values of the typical velocities. 

A similar conclusion, that the BLR radius is larger in NLSyls, was reached by Wandel 
(1997) who assumed that the representative radius of the BLR is determined by the 
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requirement to have a standard value of the ionization parameter. He then showed that 
the size of the BLR region is dependent not only on the luminosity of the central source, 
but also on the soft X-ray spectral slope. A steeper (softer) X-ray spectrum has a stronger 
ionizing power and hence, for a constant ionizing parameter, the BLR clouds are at larger 
distances from the central source, have smaller velocity dispersions and as a result form 
narrower emission lines. Laor et al. (1995) also reach a similar conclusion but in their 
picture the narrow lines in NLSyl result purely from the lower black hole mass. In our 
scenario the lower black hole mass and the shape of the SED (i.e. the steeper soft-X-rays 
which decrease T h ) combine to produce the narrow lines. 

Combining equations (l)-(4) we estimate the cloud density: 

L T 7/2 

Pc-^x^ (6) 

IVI bh ± c 

or using logarithms: 



log p c ~ log L - 2 log M bh + 7/2 log T h - log T c (7) 

As has been argued in Section 3.2.1, NLSyls have bolometric luminosities comparable 
to QSOs, although their central black holes have lower masses. The median value of a 
black hole in quasars is ~ 1O 1O M (see Czerny, Witt & Zycki 1997 calculations), while 
the median black hole mass in NLSyls (as inferred from our calculations, using the same 
ADC model - see Table 6) is 1O 8 - 26 M i.e. ~ 55 times lower. Let us assume that a typical 
quasar SED is composed of a power law and an accretion disk spectrum peaking at 10 eV 
(logi/ = 15.38, 1240A), while a typical NLSyl SED has a power law and a disk peaking 
at 80 eV (however note that the most extreme NLSyl RE J1034+396 has its peak at 
120eV - see Puchnarewicz et al. 1995). Krolik and Kallman (1988) calculated the Compton 
temperatures of the hot phase for these SEDs, normalizing both to have the same total 
ionizing energy. The lOeV bump spectrum gave Compton temperatures ~ 3.0 xl0 7 K while 
the 80eV bump gave a lower temperature ~ 8.0 xl0 6 i^. At the same time the temperature 
of the cool phase increased by a factor ~ 3.0 (0.5 in logarithm see Krolik & Kallman 1988 
Fig. 2). Hence the Compton temperature of the hot phase in NLSyl and QSOs differs 
by: logT hjNLSy i — \ogT ht Q SO = —0.57 and the temperature of the cold phase is larger by: 
logT Cj jVLs y i — logT Cj Q5o = 0.5. Substituting the above values into equation (7) implies that 
\°SPc,NLSyi ~ l°g Pc,qso ~ 1 i- e the densities of the BLR should be higher by a factor of 10 
in NLSyl than in typical QSOs with redder BBB. 

The larger BLR radii and larger by a factor 10 densities obtained from our modeling 
(as being due to hotter BBBs) are consistent with the narrow lines and line ratios observed 
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in NLSyls. Thus we conclude that the unusually hot and strong BBB in NLSyls can 
naturally produce their observed UV spectra. 

4. NLSyl vs. BAL QSOs 

It has been suggested (e.g. Leighly et al. 1997, Lawrence et al. 1997) that there 
may exist a connection between NLSyl and BAL QSOs. Both these classes have strong 
FeIIA4570 and A1IIIA1857 emission and weak CIVA1549 and [OIII]A5007. Their continua 
are red in the optical and strong in the IR; additionally both classes are mostly radio-quiet. 
Leighly et al. (1997) reported evidence for relativistic outflows in three NLSyl. 

Observationally there are also many differences. NLSyl are strong soft-X-ray emitters, 
while BAL QSOs are weak, possibly due to X-ray absorption (Mathur, Elvis & Singh 1995). 
BAL QSOs are thought to be seen more edge-on, at viewing angles skimming the edge of 
the dusty torus (Turnshek et al. 1996, Aldcroft, Elvis & Bechtold 1993). NLSyls, on the 
other hand, are probably viewed more face-on, as they show low absorption from the torus 
(Boiler, Brandt & Fink 1996) and some even show beaming in their radio spectra (e.g. 
PKS 0558-504, Remillard et al. 1991). 

In high resolution HST spectra NLSyl show absorption features which are much 
weaker than in BAL QSOs (see Table 2). However this is expected since 50% of Seyfert Is 
show absorption features (Crenshaw et al. 1995). The optical spectra of some BAL QSOs 
may resemble the spectra of NLSyl, showing narrow FL3 with FWHM < 2000 kms -1 (that 
is why the two BAL QSOs: PG 1351+640 and PG 1411+442 were initially chosen to be in 
our sample), but this only cautions us that basing classifications on optical spectra alone is 
potentially misleading. 

5. Conclusions 

In this paper we have studied the UV emission line properties of a class of extreme 
opt/UV/X-ray AGN: the narrow line Seyfert 1 galaxies. We found 11 NLSyls that had 
been observed in the UV by either HST or IUE. We have shown that in comparison with 
"normal" broader line AGN, the equivalent widths of CIV and Mgll are significantly 
smaller (NLSyl have EW(CIV)<60 and EW(MgII)<20, normal AGN have EW(CIV)<210 
and EW(MgII)<120), the EW of A1III larger (few A), and the UV line widths are narrower 
(although not as narrow as the optical H/5 line). Also the CIIIj/Lya;, CIV/Lya and 
Mgll/Lyo; line ratios are smaller, while those of Silll]/ CIII] , SiIV+OIV]/CIV lines are 
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larger. Photoionization models predict that these line ratios are formed in material with 
densities higher, by a factor few (less than 10) than standard BLR cloud densities, and 
with the ionization parameter lower by a factor 10. These parameters however predict 
higher Mgll/Lyo; ratio, in contradiction to the lower ratios observed requiring that Mgll be 
produced in a separate region. 

We have fitted the SEDs of our NLSyls to the Witt, Czerny & Zycki (1997) model 
of an accretion disk with a Compton cooled corona and found that NLSyls radiate at 
0.27 < L/LEdd < 0.58, much larger than the typical AGN (L/LEdd < 0.3). The masses of 
the central black holes calculated from the model are, in our objects, of the order of masses 
found in typical Seyfert 1 galaxies (1O 8 M ) but the bolometric luminosities [yL v ~ 10 46 erg 
s _1 ) are larger and comparable to those of QSOs. 

Krolik & Kallman (1988) predict that steeper soft-X-ray BBBs, such as these of 
NLSyls, change the equilibrium of the two-phase cloud-intercloud medium, decreasing the 
temperature of the hot intercloud medium (which we assume to be the corona above the 
accretion disk) and increasing the temperature of the cool BLR clouds. We show that this 
change in equilibrium increases the density of the BLR clouds resulting in a change of 
the observed line intensities and ratios consistent with these in NLSyls. In addition the 
resulting decrease in Th, causes an increase in the radius of the BLR, a correspondingly 
lower velocity dispersion and narrower lines as observed in NLSyls. 

The NLSyls lie at the extreme end of the Boroson and Green eigenvector 1 (Boroson 
& Green 1992), which was then found (Brandt & Boiler 1998) to link the soft X-ray 
properties with the optical properties i.e. the FeII/H/5 and [OIII] strengths and Hf3 line 
width. We have found that the NLSyls have very weak CIV and CIII] lines, and narrow 
UV lines extending the set of parameters linked to eigenvector 1. The large BLR cloud 
densities, deduced from these characteristic UV line ratios, are probably due to the steep 
soft X-ray SEDs, which are in turn, the result of larger L/LEdd ratios (as inferred from the 
Witt, Czerny & Zycki 1997 ADC model). In this scenario a larger L/L E dd is the physical 
parameter driving the Boroson & Green eigenvector 1. 

We are grateful to Niel Brandt for helping us to obtain a complete list of known 
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discussions and thank the anonymous referee for comments that improved the manuscript. 
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Table 1. Sample 



Name 


a(J2000) 


<5(J2000) 


z 


log vL v 
2500A a 


IZwl* 


00 53 34.94 


+12 41 36.2 


0.0611 


44.89 


E 0132-411 


01 34 57.36 


-40 56 22.4 


0.266 


45.00 


NAB 0205+024 


02 07 49.86 


+02 42 55.9 


0.1564 


45.24 


PKS 0558-504 


05 59 47.37 


-50 26 51.8 


0.137 


45.31 


PG 1211+143* 


12 14 17.60 


+14 03 12.5 


0.085 


45.06 


IRAS 13349+2438 


13 37 18.73 


+24 23 03.3 


0.107 


45.11 


PG 1351+640 BAL 


13 53 15.78 


+63 45 44.8 


0.087 


44.99 


PG 1411+442* BAL 


14 13 48.39 


+44 00 13.6 


0.090 


44.96 


MRK 478* 


14 42 07.46 


+35 26 22.9 


0.0781 


44.85 


PG 1444+407* 


14 46 45.95 


+40 35 06.0 


0.267 


45.44 


IIZwl36 


21 32 27.81 


+10 08 19.5 


0.061 


45.04 



* HST spectrum 

BAL: also classified as a broad-absorption line BAL QSO 

a: the luminosity at 2500A has been obtained by extrapolating B (or V) magnitude and 
assuming a continuum slope of 0.5 (H = 50 km s _1 Mpc -1 ) 
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TABLE 2 
UV Line Rest Frame EW in A" 









Emission lines 








Absorption 


lines 


Name 


Lya 


CIV 


SilV+OIV] 


Silll] 


cm] 


Aim 


Mgll 


Lya 


NV 


CIV 


IZwl* 


171.4 


27.9 


17.2 


9.8 


3.3 


7.6 


25.0 




1.4 




0132-411 6 


<140.8 


<49.5 


















NAB 0205+024 


60.9 


38.5 














0.9 




PKS 0558-504 


32.9 


13.7 


<5.4 
















PG 1211+143* 


160.8 


49.2 


8.4 


3.5 


4.9 


0.0 


18.4 




1.0 




IRAS 13349+2438 6 


<85.9 




















PG 1351+640 BAi 


140.2 


>30.0 


25.5 








12.9 




< 5.2 


>2.1 


PG 1411+442* B - 4L 


>79.5 


43.6 


12.4 


5.4 


15.8 


2.0 


18.7 


>1.0 


14.0 


7.5 


MRK 478* 


123.0 C 


36.0 


19.5 C 


9.9 


6.6 


4.7 










PG 1444+407* 


125.5 C 


25.8 


22.7 


4.6 


7.5 


4.1 










IIZwl36 


110.4 


56.6 


23.3 








32.0 









* denotes HST spectra 

BAL: also classified as a BAL QSO 

a: typical error ~ 15% 

b: low S/N spectrum 

c: equivalent width obtained from IUE spectrum 
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TABLE 3 
Like ratios 



Name 


CIV/Lya 


CIII]/Lya 


CIII]/CIV 


SiIII]/CIII] 


SiIV+OIV]/CIV 


SiIV+OIV]/Lyo 


Mgll/Lya 


IZwl* 


0.11 


0.01 


0.11 


2.97 


0.85 


0.09 


0.07 


0132-411 


>0.21 














NAB 0205+024 


0.41 


- 


- 


- 


- 


- 


- 


PKS 0558-504 


0.31 


- 


- 


- 


<0.42 


<0.13 


- 


PG 1211+143* 


0.29 


0.03 


0.08 


0.71 


0.18 


0.05 


0.03 


IRAS 13349+2438 
















PG 1351+640 Bil 


>0.16 








0.26 


0.04 


0.07 


PG 1411+442* BiI 


<0.30 


<0.13 


0.35 


0.34 


0.30 


<0.09 


<0.07 


MRK 478* 


0.24 


0.04 


0.15 


1.5 


0.57 


0.14 




PG 1444+407* 


0.17 


0.04 


0.22 


0.61 


0.87 


0.15 




IIZwl36 


0.34 








0.43 


0.15 


0.01 


Observed range 


0.11-0.41 


O.OK0.13 


0.08-0.35 


0.34-2.97 


0.18-0.85 


0.04-0.15 


0.01-0.07 


Mean 


0.25+0.09° 


0.05+0.05 


0.18+0.11 


1.23+1.07 


0.49+0.26 


0.11+0.04 


0.05+0.03 


Seyfert Is obs. range 1 


0.35-2.01 


0.03-0.39 


0.08-0.5 








0.07-0.63 


Seyfert Is mean 1 


0.53 




0.20 








0.17 


QSOs observed range 


0.3-1.04= 


0.15-0.3= 


0.46 d 


0.3+0.1" 


0.3 d 


0.08-0.24= 


0.15-0.35= 


LOC maximum reprocessing' 


0.54 


0.28 s 






0.15 


0.08 


0.38 


LOC inegrated' 


0.57 


0.12 s 






0.11 


0.06 


0.34 



* HST spectrum; BAL - also classified as a BAL QSO 
a: standard deviation a = ^ ^".jf 1 * ~ *) 2 
b: Wu et al.(1983) 

c: Baldwin et al. (1995) and Wilkes et al. (1998) 
d: composite spectrum from Francis et al. (1991) 
e: Laor et al. (1995) 

f: line ratios from Baldwin et al. (1995): LOC maximum reprocessing - calculated considering emission from clouds with optimal density and ionizing flux parameters for each line; LOC 
integrated - calculated integrating over clouds with all posible density and ionizing flux parameters and using the distribution of clouds on the density-flux plane (the assumed dependence 
with radius and density was: /(r) oc const and g(r) oc n _1 respectively) as a weighting function. 

g: CIII]+SiIII]+AlIIII/Lya 
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Table 4. Line widths in km/s 



Name 




/"ITT T 

UlV 


TV T „TT 

Mgll 


T T Qn 

tip 


IZwl* 


1730 


3190 


1850 


1240 


0132-411 a 


<3920 


<2460 




1930 


NAB 0205+024 


1630 


2520 




1100 


PKS 0558-504 


3820 


3540 




1500 


PG 1211+143* 


1600 


1940 


1945 


1900 


IRAS 13349+2438 6 


<7000 






2100 


PG 1351+640 BAL 


2370 


>2110 


2140 


860 


PG 1411+442* BAL 


>1850 


2300 


1780 


2670 


MRK 478* 


2810 


2820 




1370 


PG 1444+407* 


3700 


6180 




2480 


IIZW136 


2130 


2363 


1980 


2060 



Mean 2960+1611 2942+1237 1939+133 1746+577 



QSOs c 5399+2757 4793+1765 5566+2426 

QSOs d - 5150+1680 4580+1890 

QSOs e 3454+1291 4335+1550 3774+2212 



* HST spectra; BAL - also classified as a BAL QSO 

a: after Boiler, Brandt, Fink (1996) 

b: very low S/N spectra 

c: Wilkes et al. (1998) 

d: Baldwin, Wampler & Gaskell (1986) 

e: Corbin & Boroson (1996) 
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Table 5. Observed soft and hard X-ray indices 



Name 


a(0.1 - 2.5keV) 


Ref. 


a(2 - lOkeV) 


Ref. 


IZwl 


2.0+0.1 


1 


1.3+0.06 


7 


0132-411 


3.1+0.6 


1 






NAB 0205+024 


2.8+0.5 


2 


1.09+0.10 


10 


PKS 0558-504 


1.9+0.1 


3 


1.26+0.05 


3 


PG 1211+143 


2.13+0.22 


4 


1.01+0.06 


8 


IRAS 13349+2438 


1.70+0.33 


5 


i 99+O.O8 
± - zz -0.08 


8 


PG 1351+640 BAL 


1.53+0.63 


4 






PG 1411+442 BAL 


1.97+0.46 


4 






MRK 478 


2.6+0.1 


1 


0.95+0.13 


8 


PG 1444+407 


2.2+0.3 


1 






IIZwl36 


1.25+.. . 


6 


1.17±0.07 a 


9 



BAL - also classified as a BAL QSO 

a: slope from Ginga determined between 2-18keV (ref. 9) 

References: 

1 - Boiler, Brandt k Fink (1996), 

2 - Fiore et al. (1995), 

3 - Brandt, private communication, 

4 - Wang, Brinkmann & Bergeron (1996), 

5 - Brandt, Fabian & Pounds (1996), 

6 - Wang, Lu & Zhou (1998) (error on a(0.1 — 2.§keV) not available), 

7 - Hayashida (1997), 

8 - Brandt, Mathur & Elvis (1997), 

9 - Lawson & Turner (1997), 

10 - Fiore et al. (1998) 
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Table 6. Soft and hard X-ray indices from best fit model 



Name 


L/ L E dd 




log My, 


a< softX 


a hardX 


IZwl 6 


0.58 


0.027 


7.93 


2.02 


1.17 


0132-411° 


0.27 


0.30 


8.29 


2.52 


1.05 


NAB 0205+024 


0.27 


0.30 


8.45 


2.50 


1.04 


PKS 0558-504 6 


0.58 


0.027 


8.89 


1.89 


1.14 


PG 1211+143 


0.30 


0.14 


8.26 


2.02 


1.05 


IRAS 13349+2438 


0.58 


0.03 


8.12 


1.99 


1.15 


PG 1351+640 BAL - a 


0.27 


0.03 


8.20 


1.51 


1.11 


PG 1411+442 BAL - a 


0.58 


0.027 


7.96 


2.00 


1.17 


MRK 478 


0.27 


0.3 


8.12 


2.53 


1.07 


PG 1444+407 a 


0.27 


0.3 


8.61 


2.49 


1.04 


IIZwl36 


0.27 


0.03 


8.20 


1.51 


1.11 



BAL - also classified as a BAL QSO 

a: these objects do not have observed hard-X-ray slopes, hence the derived model 
parameters are not well constrained 

b: we could not fit both the soft and hard X-ray slopes simultaneously 

c: soft X-ray index measured from 0.1-2.5 keV corresponding to a ROSAT slope 

d: hard X-ray index measured from 2-10 keV corresponding an ASCA slope 



